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Abstract 


A new method is presented for obtaining indefinite integrals of com¬ 
mon special functions. The approach is based on a Lagrangian formnla- 
tion of the general homogeneous linear ordinary differential equation of 
second order. A general integral is derived which involves an arbitrary 
function, and therefore yields an infinite number of indefinite integrals for 
any special function which obeys such a differential equation. Techniques 
are presented to obtain the more interesting integrals generated by such 
an approach, and many integrals, both previously known and completely 
new are derived using the method. Sample results are given for Bessel 
functions. Airy functions, Legendre functions and hypergeometric func¬ 
tions. More extensive results are given for the complete elliptic integrals 
of the first and second kinds. Integrals can be derived which combine 
common special functions as separate factors. 

1 Introduction 

In tables of integrals such as [1-4], the bulk of the content consists of definite 
integrals, with almost token collections of indefinite integrals. However, the 
relatively few indefinite integrals which have been published tend to be well 
known and heavily used. This paper presents a new and surprisingly simple 
method for deriving indefinite integrals of both elementary and special functions, 
provided the function satisfies an ordinary linear differential equation of the form 


y" {x) + p (x) y' (x) +q(x)y (x) = 0. 


( 1 ) 


The principal result derived here is the indefinite integral 
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( 2 ) 


f{x) {h' {x)y{x) -h{x) y' {x)) 
where y (x) is any solution of equation o, / (x) is given by 

/ (x) = exp p (x) dx^ (3) 

and h (x) is an arbitrary function. As h (x) is arbitrary, this equation yields an 
infinite number of integrals for any function y {x) which satisfies equation (HD. 
Many well-known integrals in the literature can be derived very simply from 
equation Q, together with a large number of interesting new integrals. Tech¬ 
niques for exploiting equation ([5]) to obtain interesting integrals are presented 
here, together with sample results for selected special functions. Somewhat 
more than sample results are presented for the complete elliptic integrals of the 
first and second kinds, but it is impossible to do more than scratch the surface 
for any particular function in a single paper. It seems that the total number of 
indefinite integrals for special functions can be increased by a large factor using 
equation ([2|). In particular, it is possible to derive integrals combining products 
of different special functions, such as Bessel and Legendre functions combined 
with each other, or with elliptic integrals, or other special functions. Equation 
m was originally derived from an Euler-Lagrange equation, but once known, it 
can be proved in an elementary manner without variational calculus and both 
derivations are given here. 

Section 2 below presents the Euler-Lagrange formulation of equation (HD and 
the derivation of equation ([^D from it, together with a simpler proof. A second 
integral is also derived, which will not be considered in detail here. Section 3 
presents techniques for exploiting equation ([2]) and some well-known results are 
derived using this method. It is shown that functions which are conjugate, in the 
sense that the differential equations they satisfy have the same p (x) in equation 
HD, can always be combined as products in the same integral. It is shown that 
any two equations can be transformed to be mutually conjugate and transfor¬ 
mation equations are given to calculate q (x) for arbitrarily specified p (x). It is 
shown that the reverse process, finding a suitable p (x) to give a desired q (x ), is 
governed by a Riccati equation [5]. Both known and new integrals are derived 
in this section as examples of the methods presented. Section 4 presents sample 
results for a selection of special functions. Section 5 gives somewhat more than 
sample results for the complete elliptic integrals of the first and second kinds, 
using the method of fragmentary equations presented in Section 3. Table 1 gives 
the special functions used. 


2 Formulation 


Equation HD can be expressed in Lagrangian form as: 


dx \dy') 
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(4) 


Symbol 

Special Function 

Ai (x) 

Airy function of the first kind 

Bi (x) 

Airy function of the second kind 

Ai' (x), Bi' (x) 

Derivatives of the Airy functions 

Eik) 

Complete elliptic integral of the second kind 

2 F 1 (a, 6; c; x) 

Gauss hypergeometric function 

Gi (x), Hi (x) 

Scorer functions 

In ix) 

Modified Bessel function of the first kind 

J 1/ ix) 

Bessel function of the first kind 

Hn ix) 

Struve function of the first kind 

Kik) 

Complete elliptic integral of the first kind 

K„ ix) 

Modified Bessel function of the second kind 

Pv ix) 

Legendre function of the first kind 

PIf (x) 

Associated Legendre function of the first kind 

Qu ix) 

Legendre function of the second kind 

Qi) (x) 

Associated Legendre function of the second kind 

ix) 

Lommel function 

K ix) 

Bessel function of the second kind 

r(x) 

Gamma function 


Table 1: Special Functions Used 


where the Lagrangian L is given by 

L = fix) (x) - q (x) (x)) 

and the function / (x) is a solution of the equation: 

f'(^) 


f(x) 


= Pix) 


(5) 


( 6 ) 


and hence is given by the integral ([3]). A principal property of any Lagrangian 
is: 

d 


— [ L- ^ = — 
da; V ^ dy' j dx 


and from equation ([S]) we obtain 
,dL 


L - y'^ = -f (x) {y'^ ix) + q (a:) (a;)) 


and 


dx 


= f (a:) y' (x) - if (x) q ix)) y ix). 


(7) 


( 8 ) 


(9) 


Substituting equations dH) and dH) into equation dZ]) gives immediately the in¬ 
tegral 

[ {if (a;) q (a;))' y^ (a;) - /' (a;) y'^ (a;)) dx = f (x) {y'^ (x) +qix)y ix)) (10) 
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which is analogous to the energy integral in classical mechanics. This integral 
can be considered to be complementary to equation as it gives different inte¬ 
grals, but these results will be presented separately. Explicitly evaluating both 
derivatives in equation © using equations and m gives after cancellation 
and collection of terms: 

-if (x) y'(x)) = f {x) q (x) y (x) (11) 

and we obtain a second indefinite integral: 

J f (x) q (x) y (x) dx = -f (x) y' (x). (12) 


The differential equation ([5]) can be transformed by defining y{x) = h (x) z (x) 
to give: 

h" (x) h' {x) 




h ((x)) h (x) 


p(x) + q(x) j z = 0. (13) 


This transformed equation can also be put in the Lagrangian form 


dL 


dL 


where 

and 


dx I dz' ) dz ^ 


L = f{x) (z'2 (x) -qix)z (x)) 

_ h" (x) h'{x) 

d i.x) = , , ■! + -tTTP (^) +d{x). 


h (x) h (x) 

The function / (x), the equivalent of / (x) in equation ([5]) , is given by 

h' (x) 


/ (x) = exp 


h (x) 


-I- p (x) ) dx I = (x) / (x) 


so that the new Lagrangian can be expressed as 


L = h‘^ (x) / (x) ( - 


Pix)+q (x) z^ (x) . 


h" (x) h' (x) 
h{x) h (x) 

The equivalent of the integral m for the transformed equation is 


(14) 

(15) 

(16) 

(17) 

(18) 


J (x) / (x) (^) + d 2 (x) dx 

and since z (x) = y (x) /h (x) and hence 

, y'{x) h'{x)y{x) 


(x) / (x) z' (x) 

(19) 

( 20 ) 


then equation (HI ) reduces to equation @. This relation holds for all ho¬ 
mogeneous second-order linear ordinary differential equations for an arbitrary 
function h (x). Once known, this result can be proven very simply without 
reference to variational calculus. 
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Theorem 1 


f{x){h'{x)y{x)-h{x)y'{x)) (21) 

where f (x) and y (x) obey the respeetive differential equations: 

f {x)=p{x)f{x) (22) 

y" {x)+p{x)y'{x)+q{x)y{x) =Q (23) 

and h(x), p (x) and q {x) are arbitrary funetions. 

Proof. Arbitrary functions / (cc), h (x) and y (x) satisfy the differential identity: 
^ (/ (a;) (h' (x) y{x)-h {x) y' (x))) = /' (x) {h' {x) y{x)-h (x) y' (a;)) (24) 

+/ {x) {h" (x) y{x)-h {x) y" (x)). 

Eliminating /' (x) and y" (x) from equation ([24|) using equations (l2^ and (l23ll . 
respectively, gives 


^ (/ (x) {h' (x) y{x)-h {x) y' (x))) = 

/ (x) {h" (x) + h' (x) p{x) + h (x) q{x))y (x) (25) 

for arbitrary h (x) and integration of equation (l25l) gives and the theorem 
is proven. Equation (l22l) can be integrated immediately to give 

f (x) = exp p {x) dx'^ . (26) 


In equation ([2]) we can take y (x) = yi (x) to be any solution of equation 
m and h (x) = j /2 (x) to be another solution, to obtain the Wronskian for the 
differential equation o as: 


W(x) 


yi (x) 2/2 [x) 
y'l (a;) 2/2 {x) 


A 

fix) 


(27) 


where A is a constant. Erom equation ([3]) this is equivalent to Abel’s formula 
[6,7] for the Wronskian: 


W {x) = W (xq) exp 





(28) 
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3 Techniques for finding indefinite integrals 

Equation Q yields an infinite number of indefinite integrals for each solution 
y (x) of equation 0, as any h (x) yields an integral. In equation Q we can 
take h (x) = 1 or substitute any elementary function such as e“^ or and 

obtain an integral for any special function. Expressions involving elementary 
functions, such as h (x) = x"* and h (x) = often yield integrals which 

are interesting and new. We can also specify h (x) to be the solution to a 
differential equation where one or two terms in equation Q is deleted, such as: 

h" (x) + p (x) h' (x) = 0 (29) 

p (x) h' (x) q{x)h (x) = 0 (30) 

h" {x) + q{x)h{x) =Q. (31) 

It is frequently the case that p (x) or q (x) may consist of several terms, and 
h (x) can be specified as the solution of a differential equation where some of 
these terms have been deleted. These equations with one or more terms deleted 
will be referred to as fragmentary equations. 


3.1 Integration of the fnnction y (x) itself. 

To integrate a solution y (x) of equation ((2) it is necessary to specify h (x) to 
be a solution of the inhomogeneous equation: 


h" (x) + p (x) h' (x) q{x)h (x) = 


1 


fix) 


For the Legendre equation: 


2x 


then 

/ (x) = exp 1^- J 
and equation (15^ becomes 


2x 


1 — x' 


7 // / N 2x ,, , s ly (i^ + 1), , . 1 

^ “ 1- 2 ^ + “i-2^^ = 1- 

1 — x^ 1 — x'^ 1 — c 


which has a simple solution 


h (x) = 


iy{iy + l)' 

Substituting (1551) into equation ([5]) for Legendre’s equation gives 


dx = 


x^ — 1 


jPLix) 

Qi7(x)J“'“' iy{iy+ 1 )\Q'^{x) 


Pi' (x) 


(32) 


.<‘'+)bw = o 

1 — 

(33) 

-Trdx ) = 1 — X^ 

7 / 

(34) 


(35) 


(36) 


(37) 
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and the Legendre recurrence relations [1]: 



\q„+i (x) - xQ^ (x) J 

_ ^ f XP^ (x) - Py-i (x) 1 
\xQi, (x) - Q^-i (x)/ 


then give the alternative forms 


/■ I Pj. (x) 1 _ 1 f Pi,+i (x) - xP^ (x) 1 

J \Q^{x)j i^\q„+i (x) - xQ„ (x)) 


(38) 

(39) 


(40) 


/■ f Pi, (x) _ 1 [ xP^ (x) - Pi,_i (x) 1 


(41) 


The Legendre equation above is a simple case, and sometimes the solution to 
equation (I32p is given in terms of functions specifically defined to satisfy this 
inhomogeneous equation . Examples are the Lommel [8,9,] and Struve [9] func¬ 
tions for Bessel’s equation and the Scorer functions for Airy’s equation [2,10]. 
The Bessel equation is 


y" (x) -h -y' (x) + (l - y {x) = 0 
X \ J 


(42) 


and has / (x) given by 

/(x) = exp = X. (43) 

The Lommel function Sm,n (x) satisfies the equation 

y" (x) + -y' (x) (l - y (x) = x™“^ (44) 

X \ x^ J 

with equation (15^ being the special case for m = 0. Taking h (x) = Sm,n (x) in 
equation @ gives 

j X™ J„ (x) dx = X (s(„ „ (x) Jn (x) - Sm,n (x) J'^ (x)) (45) 

which from the Lommel and Bessel recurrence relations [1]: 

Tl 

n (^) — Tl l) Sm—l,n—l (^) Jn (^) ^m^n (^) Jn (^) 

’ X 

Tl 

-Jn (x) + jy (x) = Jn-1 (x) (47) 

X 

reduces to 

X {s'^,n (x) Jn (x) - Sm,n (x) (x)) (48) 
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and gives 


/ X^Jn (x) dx = X ((to + n - 1) Sm-l,n-l (x) Jn (x) - Sjn,n (x) J„-l (x)) . 

(49) 

Watson [9] gives a rather different derivation of this result. For the special case 
where m = n, substituting the identities [1] 

Sn,u {x) = r (n + 1/2) (x) (50) 

into equation (1491) gives the integral in terms of the Struve function H„ (x) as: 


f x^Jn (x) dx = v^2” (n + 1/2) x (H„_i (x) J„ (x) - H„ (x) J„_i (x)). 

(52) 

The Airy functions are a similar case to that of the Bessel functions described 
above. The Airy functions Ai (x) and Bi (x) obey the differential equation: 

y" (x) - xy (x) = 0 (53) 


and the Scorer functions Gi (x) and Hi (x) obey the inhomogeneous Airy equa¬ 
tions [10]: 

y" (x) - xy (x) = T- (54) 

TT 

where the negative sign is taken for Gi (x) and the positive sign for Hi (x). 
Taking h (x) in equation @ to be either Gi (x) or Hi (x) gives the alternative 
integrals, which of course differ only by a constant: 


'Ai (x)^ 
,Bi (x). 

^dx = Tr< 

(Ai' (x) Gi (x) — Ai (x) Gi' (x) 

[ Bi' (x) Gi (x) — Bi (x) Gi' (x) \ 

Ai (x)j 
Bi(x)J 

>dx = 7r| 

' Ai (x) Hi' (x) — Ai' (x) Hi (x) 1 
^Bi (x) Hi' (x) — Bi' (x) Hi (x) / 


(55) 


(56) 


3.2 Conjugate differential equations 

Two differential equations of the form 0 are considered conjugate if they have 
the same p (x), and hence the same / (x), but different q (x). For example, the 
equation 

y" (x) +-y'(x) + ^ y{x)=0 (57) 

X 1 — x^ 

is conjugate to the Bessel equation and has as a general solution [1] the 
complete elliptic integral expression 

y (x) = CiE (x) + C 2 (E (x') - K (x')) (58) 
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where x' = \/\ — is the complementary modulus. The Bessel equation 
can be slightly generalized as 

y” (x)+ -y'(x) + (a^y{x) =0 (59) 

X \ x^ J 

which has the solution y (x) = CiJn {ax) + C 2 Yn {ax) and two equations of the 
form ()59ll are mutually conjugate unless both a and n are the same for both. 
The modified Bessel equation 

y" {x)+ -y'{x) - (a^+ ^^^y{x) =0 (60) 

X \ x^) 

with solution y {x) = Ciln {ax) + C 2 Kn {ax) is always conjugate to equation 

(IMl). 

In equation ([2]) we can take h {x) to be any solution of a conjugate equation 
with q {x) instead of q {x). We have 

h” {x) + p {x) h' {x) = —q {x) h (61) 

and hence 

h" {x) + p {x) h' {x) + q{x) h {x) = {q {x) — q{x))h {x) (62) 

so that equation ([2]) gives the integral 

J f {x) {q {x) - q (x)) h (x) y (x) dx = f (x) {h' (x) y {x) - h (x) y' (x)) (63) 


If in this equation we take y (x) = (ax), where (ax) = aJ„ (ax) + (x) 

is any solution of the Bessel equation (l5^ , and h (x) = {l3x) is any solution 

of the Bessel Equation 


h" (x) + —h' (x) + f ^ ) h (x) = 0 


X \ x^ 

then equation (IS51) gives the well-known Bessel integral [1]: 

y ^(a^ - /3^) X - -— ^ Zn {ax) Zm (/3x) dx = 


(64) 


X (/3^m iPx) Zn {ax) - Zm (/3x) Z'^ (ax)). (65) 


Similarly, taking y (x) to be any solution of equation (1591) and h (x) = Zm {l3x ), 
where Zm (/3x) = Ain (/3x) + BKn (/3x) is any solution of the modified Bessel 
equation: 


h" 


(x) -I- -h' (x) -I- 

X 



h{x) = 0 


( 66 ) 
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gives the integral 


2 


/ 



Zn [ax) Zm [Px) dx = 


X {pz'^ [Px) Zn [ax) - Zm [Px) Z'n [ax)) . (67) 

In equation ((63ll we can take y [x) = (ax) and (x) = E (x), a solution of 
the conjugate equation ((57ll . to obtain 




dE (x) 


Zn [ax) — E (x) 


Zn [ax) E (x) dx = 
dZn [ax) 


where [1]: 


dx ' (^ 2 ; 

dE (x) _ E (x) - K (x) 
dx X 

A particularly simple special case of equation (l69ll is 


( 68 ) 

(69) 


Vo [x) E (x) dx 


Jo (x) (K (x) - E (x)) - X Ji (x) E (x). 


(70) 


3.3 Transformations of the differential equations 

Any differential equation of the form o can be made conjugate to any other by 
a change of dependent variable. This allows integrals to be obtained containing 
a product of any two specified special functions, provided they both satisfy 
differential equations of the form ©■ The differential equation 

z" (x) + p (x) z (x) + q[x) z (x) = 0 (71) 

can be made conjugate to equation m by the transformation y[x)=g (x) z (x) 
which yields the equation 


y" [x)+ 


2g' (x) 


9[x) 

and we can choose 

so that 


+ p(x) ) y' (x) + 

[x) 
a[x) 


g" [x) ^ g' [x) 
g [x) g [x) 

+ p[x) =p[x) 


P [x) +q[x)]y (x) 


g (x) = exp Q J [p (x) - p (x)) dx^ . 
From equation (TM)) we obtain the relation: 


g[x) = 


' f_[x) 

f[x)' 


0 (72) 

(73) 

(74) 


(75) 


10 













The coefficient of y (x) in equation dZl]) can be evaluated directly from g (x) but 
it is usually more convenient to note that 


d (g'{x)\ g"{x) (g'{x)V 

dx\g{x)j g (x) [ g (x) J 

and 

^-7^ = ^ (P (x) - P (x)) (77) 

g(x) 2 

so that equation dzii) becomes 

y" (x) + p (x) y' (x) + (p (x) - p (x))' + P ^ P (^) q y ( 3 ^) ^ q. 

(78) 

For example, the equation satisfied by the complete elliptic integral of the first 
kind K {x) is 

which has f (x) = a: (l — x^) and general solution [1]: 

z{x) =CiK{x) + C2K{x'). (80) 


The transformed equation conjugate with the Bessel equation is given by equa¬ 
tion ([7^ as 

y” (x) + -y' {x) + ^ ^ y {x) = 0 (81) 

X (1 — X‘‘) 

and from the relation (175]) the general solution of equation (|5T|) is given by 

y (x) = Cix'K (x) -b Cax'K (x'). (82) 


In equation (IS51) we can take y (x) to be (ax), the general solution of the 

Bessel equation (l5^ and specify h (x) = x'K (x). This gives the integral 


/ 


xx' 




Zn(ax) K(x) dx = 


X (^Zn (ax) ^ (x'K (x)) - 

where [1]: 

dK(x) _ E(x) 
dx xx'^ 

A special case of equation (1551) is 



(83) 

K(x) 

(84) 


r x^ (2 — x^) 

J x'^ 


f Jo (x) 
\Yo (x) 


K (x) dx 


f Jo (x) 
IVo (x) 


K(x)-E(x) ,rJi(x) 
x' '''' iFi (x) 


K(x). 

(85) 
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Employing the conjugate equations (1571) and (ISTl) in equation (1551) . with y (x) = 
E (x) and h {x) = x'K {x) gives the integral 

J ^ E {x) K (x) dx = (K (x) — E (x)) ^^ — x'K (x)^ . (86) 

If the common p (x) of two conjugate differential equations is transformed such 
that p (x) —>■ p (x) and / (x) —>■ / (x) for both equations, the integral given by 
equation (1631) remains unchanged. This can be shown by applying the transfor¬ 
mations y (x) = g (x) y (x) and h(x) = g (x) h (x) to the conjugate equations: 

y" {x)+p[x)y'{x) + q{x)y{x) = Q (87) 

h" (x) -I- p (x) h' (x) -I- q (x) h (x) = 0 (88) 

to obtain the still conjugate differential equations: 

y" {x) + P {x) y' (x) -f qi (x) p (x) = 0 (89) 

h" (x) + p (x) h' (x) -I- qi (x) h{x) = 0 (90) 

where from equation (l78)) 

gi (x) = ^ (p (x) - p (x))' + ^ + g (a;) (91) 

gi = ^{pix)-p (x))' + ^ ^ +q (x) (92) 

and hence 

qi{x) - qi{x) = q{x) - q{x). (93) 

From equation (1751) the general solutions of the differential equations (l89)) and 
dSOD are related to the general solutions of equations (1571) and (1551) . respectively, 
by: 

y (a;) = \lj^y 

h (a;) = (x) (95) 

Substituting equations (l93l) - (l95l) into equation (l63l) yields after cancellation 

y/(a;) iq{x) - q{x))h{x)y{x)dx = /(x) {h' (x) p (x) -h{x)y' [x)) (96) 

and hence the transformed equations do not yield a new integral in this manner, 
though new integrals can be obtained from the transformed equations using 
fragmentary equations. 
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3.4 A Riccati Equation 

Instead of transforming p (x) in equation o, we can instead target q (x) for 
simplification by a dependent variable change. The differential equation 

z" (x) + p (x) z' (x) + q (x) z (x) = 0 (97) 

is transformed by the variable change y (x) = g (x) z (x) to give 

y" (x) + p {x) y'(x) + q (x) y (x) = 0 (98) 

and then from equation (1781) the function q (x) is given by 

q{x) = ^ {p (x) - P (x))' + ^ ^ ^ + g (a:). (99) 

If p (x) and q (x) are considered known and q (x) is to be specified, solving for 

p (x) as an unknown gives a Riccati equation [5]. Although nonlinear, it is a well- 
known equation with many useful exact solutions [5]. This can be illustrated 
with the associated Legendre equation, which has a relatively cumbersome q (x ), 
in the context of applying equation ©: 


ff { \ 

^ 1-2 ^ 

1 — 


\ I- (1-a 


z (cc) = 0. 


( 100 ) 


Under a general dependent variable change, the new q (x) is given by dMl) 


as 


1 dp p^ v{v -\-l) — 1 

2 da; 4 


1 - 


{i-x^y 


which has a simple solution 


2{x-p) 

P\^)= -1-2“ 

1 — a;"' 


and the transformed associated Legendre equation is therefore 

y “ “i ——y + “i ——y ^ 

1 — x^ 1 — x^ 


( 101 ) 


( 102 ) 


(103) 


For this differential equation / (a;) and g (x) are given by equations ® and (EH), 
respectively, as 

(104) 

(105) 




3(x) = 

and hence the general solution of (|103p is 


— X 

i-xy^^ 

1 + xj 


/ 1 — x \ 

y j (a;)) ■ (106) 
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Simple integrals can be obtained by taking h {x) to be either of the independent 
solutions of the equation 


which are 


h" ix)-‘^^^h'{x)=0 
1 — 

h{x) = 1 


h (x) = 


1 — X 
1 + X 


(107) 

(108) 
(109) 


Substituting h(x) = 1 into equation ([2]) gives 


1 + a; 

I — X 


m/2 


Pi^ (x) 

QO (x) 


dx = 


1 


1 + xy/" (x) + (x2 - 1) (a:)l 

1-xJ 1/rQO (x) + (x^ - 1) (cc) J 


iy(iy+l) yl ~ • 
which from the recurrence relation [1]: 

(1 - x^) (x) = i^xP^ (x) - (I' + pi) P^-i (x) 

can be expressed as 


( 111 ) 


1 + X 
1 — X 


m/2 


P/f (x) 
Qv (x) 


dx = 


1 f I + x\^^^ ( {vx + /r) Pjf (x) - (z/ + /r) P^_^ (x) 


z^(z^ + l) yl-x/ \(z/x + ^) QO (x) - (z/+ ^) QO-i (x) 
Substituting equation (11091) into equation ([2) gives after some reduction: 


( 112 ) 


1 — X 
1 + X 


m/2 


Plf (x) 
Qu (x) 


dx = 


1 


- xV^^ r f (z/x - /r) (x) - (z^ + /r)Pi._i (x) 


z^(z/+l) 


(irf) [{ 


_ (z/x - n) Qv (x) - (z/ + /r) Qjy-i (x) J J ^ 

Equations (I112|) and (111311 are previously known integrals which are tabulated 
in a combined form in [4]. A different integral is obtained by taking h (x) to be 
the solution of 

+ = o (114) 


1 — X^ 


which is: 


h (x) = (x —. (115) 

Substituting equation (11151) into equation ((2) gives after some reduction 

\dx = (x- f 

(1-x)^ 
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( 2 (1 - x^) fPi; (a;)1 

l^(^_l)(^ + 2) {x-fi)\Q>i {x)r 

_^_r fi)Pif{x)-{i^+n)Pi^_^{x)]\ 

{v- l)iy{v + l) {ly + 2)\{vx + fi) Q{!; (x) - {iy + n)Q'^_^ {x)jJ ' 
This integral seems to be new. 


4 Examples for some special functions 


There are an unlimited number of cases as h {x) is arbitrary. The art of using 
equation ([2]) is to choose h {x) to give an interesting integral. 

4.1 Bessel Functions 

Equation Q for Bessel’s equation gives: 

J X ^h” (x) H- h' (x) "b - 2 ^ ^ {x) dx = 


X {h' (x) J„ (a;) — h (x) (x)) 

and substituting h (x) = x™ in this equation: 

_ ...m+l ( (x) 


m+1 


1 + 


Jn (x) dx = x” 


which from the Bessel recurrence 


{x) — 2 (x) Jn+l (x)) 


(117) 


- j; (x) ( 118 ) 


(119) 


can be expressed as: 


pm+l ^ 


— n? 


Jn (x) dx = 


^m+i [^ 1 

X 2 


( 120 ) 


The second term in equation (I120p vanishes for m = ±n, and the Bessel recur¬ 
rence 

mJn (x) 1 


— 2 (^) ~b Jn+l (x)) 


then gives the elementary cases 

J X^'^^Jn (x) dx = x"^^ J „+1 (x) dx 
J X~^^^Jn (x) dx = —x“”^^ J „_1 (x) dx. 


( 121 ) 


( 122 ) 


(123) 
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The general integral given by equation (11201) appears to be new. For n = 0 and 
m = 1 this equation gives the simple expression 


y (l + x^) Jo (x) dx = xJq (x) + x^ Ji (x) 


(124) 


which can be readily verified by differentiation. 

Substituting h{x) = in equation pi7l) gives the integral 


(m^ — n^) 




sinx 

cosx 


+ {2m + 1) x" 


= X 


m+1 


mJn (x) 


- 4 (x) 


sinx 

cosx 


+ Jn (x) 


cosx 
— sinx 

cosx 


Jn (x) dx 


— sinx 


(125) 


The integrand of (jl25l) reduces to a single term for m = n or m = —1/2. 
Employing the Bessel recurrence: 


j! t'^Jn, (x) 

Jn (a^) — Jn+l (x) 


gives for m = n: 

sinx 


^cosxj 
and for m = —1/2 : 


Jn (x) 


Jn (x) dx = 


r J„ (x) sinx — J„+i (x) cosx 


{n—l/2)y/x — 1/4 


2n + 1 {Jn (x) cos X + Jn+i (x) sin x 

fx-^/^Jnix)\^^^]dx = 

J (cosxj 

^/xJn+l{x)\ J"sinx'( ^/xJn{x)f cosx 
cosx 


(126) 


(127) 


— 1/4 ( — sin 


(128) 


Equation (11271) is tabulated in [3]. Equation (11281) is also tabulated in [3] but 
unfortunately with multiple typographical errors. The same errors are also 
present in the original Russian edition [11], but it is evident that the authors of 
[3,11] did originally have equation 


Substituting h (x) = x" In (x) in equation (ITT71) gives 

J (2nx"'“^ + x"’*’^ In (x)) (x) dx = x" (1 + n Inx) Jn (x) — x^'^^ In (x) J), (x) 

(129) 

which can be simplified with equation (11261) to give 

J (2nx”“^ + x"+^ In (x)) (x) dx = x” J„ (x) + In (x) x"'*'^ J„+i (x). (130) 

For n = 0, equation (11301) reduces to 

y X In (x) Jo (x) dx = Jo (x) + x In (x) Ji (x). (131) 
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Despite its simplicity, equation (11311) appears to be new, and the Mathematica 
software [12] is unable to evaluate this integral. 

Substituting the power function h (x) = x'^ into equation @ with y (x) given 
by the Airy equation (1531) gives the integral: 


/ 


{n {n 



(132) 

Equation (11321) reduces to a single term for n = 0 and n = 1, to give the 
well-known integrals: 


Defining 




Ai {x) 
Bi {x) 


Ai(a^)'l . fAi' (x) 

Bi(cc)J iBi'(cc) 


dx = X 


Ai'(x)1 

Bi'(x)/ 


Ai (x)! 
Bi(x)j 


In ( 3 :) 


Ai (x) 
Bi (x) 


dx 


(133) 

(134) 

(135) 


then equation (11321) can be expressed as a recursion relation for the /„ (x): 


/„+3 {x) = (n + 1) (n + 2) /„ (x) - (n + 2) 

(136) 

Using equations (l5^ - (l56l) and (I133I) - (I134|) as starting formulas, equation (I136|) 
allows closed form solutions to be obtained by upward recursion for all n G Nq. 
Downward recursion to negative integer values of n is blocked when (11331) and 
(|134p are the starting formulas, but it can be done with either (|55l) or (1561) as a 
starting formula. 

The Airy functions satisfy the differential equation (1531) and from this equa¬ 
tion it follows that the slightly more general differential equation 


y" (x) - (x - a) 2 / (x) = 0 


(137) 


has the general solution 


y (x) = Cl Ai (x — a) -I- C 2 Bi (x — a). (138) 

Substituting h (x) = sin (x -I- (jd) into equation ([2]) with y (x) a solution of equa¬ 
tion (11371) for a = 1 gives 


J X sin (x -I- 4>) 


f Ai (x 
(Bi (x 



sin (x -I- (fi) 


r Ai' (x 
(Bi' (x 



cos (x -I- (/)) 


f Ai (x 
(Bi (x 



(139) 
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Similarly, substituting h {x) = into equation @ for a = —1 gives 


/ 




rAi(x+ 1) 

\Bi(a; + 1) 


dx = 


/rAi'(x + i)i rAi(x + i)n 
UBi'(cr + l)/ + lBi(a: + l)/;- 


(140) 


Substituting h (x) = Ai (x — (3) and h {x) = Ai (x — /3) into equation (I5|) gives, 
respectively, 

f fAi(x —,5)Ai(x —a)! 

J \Ai(x-^)Bi(x-a)j“^ 


1 r Ai' (x — /3) Ai (x — a) — Ai (x — /3) Ai' (x — a) 1 
a — /3 ( Ai' (x — /3) Bi (x — a) — Ai (x — /3) Bi' (x — a) / 

r rBi(x —/3) Ai (x —a)'l 
J \Bi(x-/3)Bi(x-a) 

1 f Bi' (x — /3) Ai (x — a) — Ai (x — /3) Ai' (x — a) 1 
a — /3 \ Bi' (x — /3) Bi {x — a) — Bi (x — (3) Bi' (x — a) J 

The integrals given by equations (jl39l) - (ll42l) appear to be new. 


(141) 


(142) 


4.2 Gauss hypergeometric functions 

The Gauss hypergeometric function y {x) = 2 F 1 ia,/3;j;x) satisfies the differ¬ 
ential equation [ 1 ]: 


^ , 7 - (a + /3 + 1) a: ,, ^ a/3 

y a; H----- y x)--- -y (x = 0 

x{l — X) x{l — X) 

and for this equation / (a:) is given by equation as 


(143) 


fix) 


x^ (1 - x)“+'^-^+'. 


(144) 


Substituting h (x) = 1 into equation @ gives 


J x'^ qi-x)“+^ ^ 2 ^ (a,/ 3 ; 7 ;x) 


-x^ (1 - x)“+^""'+^ 2 F 1 (a + 1, /3 + 1; 7 + 1; x) 


which is a tabulated integral [4]. 

We can take h (x) to be a solution of 


7 - (a-f l)x ,, , a/3 

■y ( 2 ^) —71—ry (®) = 0 


X (1 — x) 


X (1 — x) ‘ 


which gives 


h{x) = (■J — {a + f3 + 1) x) “+'^+1 


(145) 


(146) 


(147) 
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and substituting this h {x) into equation ([5]) gives 


(7 — (a + /3 + 1 ) a;)' 


J+Ffl' 


‘ 2 F 1 {a,/3;j;x)dx = 


x'^ (1 - 

(q;/3 + Q! + /3 + 1) 


(7 - (a + /3 + 1) x) 


_a£_ 

“ + /'+! X 


/ 2 F 1 {a,p-,^;x) _ 2-Fi(a + l,/3+l;7+l;x) 

\7 — (a + /3 + 1) X 7 


(148) 


4.2.1 Conjugate hypergeometric functions 

In equation (l2|) we can choose y (x) to be the hypergeometric function 

y{x) = 2Fi{a,(3;"f,x) (149) 

and h (x) to be a conjugate hypergeometic function: 

h{x) = 2 F 1 {a + 6, P - S-,j-,x) (150) 


where h (x) satisfies the equation 


a/3 S {a — /3 + S) 


X (1 — x) X (1 — x) 


h (x) = 0 . 

(151) 


X (1 — x) 

This gives the integral 

■ySia-lS + S) J {1 - x)°‘''~‘^~'^ 2 F 1 {a + S, ^ - 5 -,t,x) 2 F 1 {a, ^;'y;x) 


dx = 


x'*' (1 - x)""*"^ {oil32Fi {a + 5,(3 - 5] 7 ; x) 2 F 1 (a + 1, /3 + 1; 7 + 1; x) 

- (a + (5) (^ - 5) 2^1 (o + (5 + 1, /3 - 5 + 1; 7 + 1; 2 ;) 2 F 1 (a, (3-, 7 ; x)). (152) 


5 Complete elliptic integrals of the first and sec¬ 
ond kinds 

In this section the standard notation will be adopted for elliptic integrals, with 
the modulus denoted by k and the complementary modulus by k' = Vl — 
The complete elliptic integral of the first kind K (fc) is a solution of the differ¬ 
ential equation [ 1 ]: 

9k \ 1 

"" - T^) y' ^ 

which has f (k) = k (l — k^) and has the general solution 

y{k) =CiK{k) + C 2 K' (k) (154) 
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where K'(/c) = K(fc'). The corresponding equation for E(fc), the complete 
elliptic integral of the second kind, is [1]: 

{k)+ tz'{ k) + - —^z(fc)=0 (155) 

k 1 — 

which has / {k) = k and the general solution 

z (fc) = CiE (fc) + C 2 (E' {k) - K' {k)). (156) 

Here, only the solutions K (fc) and E (fc) of these equations will be considered, 
though the other solutions also lead to interesting integrals. 


5.1 Transformations of the Equations 

These equations can be transformed to give many different forms for p (fc), and 
only a few will be considered here. The dependent variable name y (fc) will be 
retained for transformations of equation (11531) and the variable name z{k) will be 
used for transformations of equation (I155|) . Equation (11531) can be transformed 
to be conjugate with both the Bessel equation (l59ll and equation (I155p . which 
gives: 

y" (fc) + iy'(fc) + (fc) = 0 (157) 

and this equation has the general solution 

y(/fc) =C'ifc'K(fc) + C2fc'K'(fc). (158) 

Transforming equation (11531) to be conjugate with the associated Legendre equa¬ 
tion (|100l) gives: 

Ok 1 

y"(k)-^,y'{k) + —y{k)=Q. (159) 

This equation has / (fc) = (l — and the general solution: 

y(fc) = C'i\/fcK(fc)-bC'2\/fcK'(fc). (160) 


Equation (11551) can be transformed to be conjugate with equation (11531) . which 
gives 

and this equation has the general solution 


z (fc) = Cl 


E(fc) 

k' 


^ E' {k) - K' {k) 

^2 - T-, - 


(162) 


Equation (11551) can also be transformed to be conjugate with equation (I159p 
and the associated Legendre equation (llOOl) . The transformed equation can be 
expressed in either of the forms: 


z"{k) 


2k 
1 - 


z' (k) + 



1 

1 - k^ 


{1-k^f 


z{k) =0 


(163) 
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«"(*)-T^"'w+( if*•<*> = “ <““> 

and has the general solution 

^ v^E(fc) , yfc(E'(fc)-K'(fc)) 

^i.k) = Ci - k’ -+- k’ -■ 

An unlimited number of similar transformations are possible, each yielding dis¬ 
tinct integrals, but only equations (I153p and (11551) will be considered in detail 
here. The technique used below is to specify h {x) as a solution to some frag¬ 
ment of the differential equation under consideration. For each equation, this 
gives a long list of integrals, most new and some very surprising. Since an un¬ 
limited number of transformations of the equations are possible, an extremely 
large number of interesting integrals can be generated. The lists of integrals 
given below for equations (11531) and (11551) are not complete, as some fragmen¬ 
tary solutions were rejected as not particularly interesting. However, the reader 
will be able to generate many additional integrals. 


5.2 Integrals from fragments of the equation for the com¬ 
plete elliptic integral of the first kind 

In equation ([21) with y (k) a solution of equation (11531) we can take h{k) = 1 
and employ the recurrence relation: 


dK (fc) _ E (fc) K (fc) 

dk ~ 1^ ~ 

to obtain the integral 

J kK (fc) dfc = E (fc) - k'^K (k). 

The six fragmentary equations: 


(166) 


(167) 



(168) 

h"{k) + yh' {k)=0 
k 

(169) 

h''ik)-Y^h'ik)=0 

(170) 

9k \ 1 

(171) 


(172) 
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—2k 1 

have the (non constant) respective solutions: 

h (k) = In ik/k') 

h (k) = In (fc) 
h (k) = arctanh (fc) 

h (k) = 


1 


(SP - 1)1/6 

1 

F 

1 

'/k 

and these give the six respective integrals: 


h{k) = 
h(k) = 


J kin K (fc) dk = In (E (k) - fc'^K (fc)) - K (fc) 


(173) 

(174) 

(175) 

(176) 

(177) 

(178) 

(179) 

(180) 


J fc (2 + In (yfc)) K (fc) dk = In (fc) E (fc) - (1 + In (fc)) K {k) (181) 

J {1 — k arctanh (k)) K (k) dk = /cK (k) — arctanh (k) ^E (k) — k' K. 


fc (1 - F) (1+ 4P) (2fc2 - l) (1 - fc2) E(fc) 

^ '-K (k) dx = 4- IA-^—Lk (fc) - ^ 


(182) 


(SP - 1)13/6 


(3P - 

fcK(/c) „ K(k)-E(k) 

-dk = - 


(3P - 1)1/6 

(183) 


A:'3 


k' 


(184) 


r k'-^ 1 

J {x) dx=^ {2k'^K (k) - 4E (fc)) . (185) 

Equations (11671) and (11841) are given in [4] but equations (I180p - (I183I) and (I185p 
appear to be new. Additional fragmentary equations can be constructed using 
the identity: 

1 e ^ ^ 

(186) 


1 - fc2 “ ^ ' 1 - fc2 
which gives equation (115311 in the alternative form: 


f,n^ /’I 2A: 


y' (k)- [1 + 




1 - fc2 


y' (k) = 0 . 


Useful fragmentary equations from equation (11871) are: 

h" (k) + i/i' (k)-h = 0 
k 


(187) 


(188) 
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h" {k)-h = Q 
1 2k 
k 


1 - A:2 

2k 


h'{k)-h = 0 


1 - fc2 
1 
k 
2k 


h'{k)-h^O 


h'{k)-h = 0 


1 - fc2 
2k 


h'(k) - 


fc2 


1 - fc2 


h = 0 


l-P 

1 


h' (fc) - 


fc2 


1 - fc2 


h = 0 




h = 0. 


These equations have the respective solutions 

hik) = CiKoik) + C2loik) 
h{k) = 

h{k) = 


(3/c2- 1)9 

7 \ C'l 1 7,2 


/i(A:) = Cie^ 




h (fc) = 




(SP- l)w 
hik) = e-i'^" 


1 

F 


e 2 


which in turn yield the respective integrals: 


2Ki (k) - kKo {k) 


—2Ii (fc) — klo (k) 

K{k)- 


K (fc) dfc = 


A:'2 (Xo (fc) - kKi (k)) 
fc'2 (Jo (fc) + fell (fc)) 


ifo(fc) 

/o(fc) 


EF) 


y (1 - 3F t e^'^K (k) dk = {{k ± 1) /c'^K (fc) T E (fc)) e 

/•A:(l-F + 6F-9fc6-fc8) 

/ —-Jg-K(fc)(ifc = 

d (3A:2-1)- 


(189) 

(190) 

(191) 

(192) 

(193) 

(194) 

(195) 

(196) 

(197) 

(198) 

(199) 

( 200 ) 
( 201 ) 

( 202 ) 

(203) 


(204) 

(205) 
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(1 - P) (fc4 + 2fc2 - 1) 


3*2 - 1 

64 


K (*) - E (*) 


ee 


( 3*2 - 1)5 


(206) 


(1 + *2 - 5*^ - *®) ^-K (*) dfc =(2(1- *'‘) K (*) - 4E (*)) 


*2 




(207) 


J * (*^ + 3*2 - 1) (fc) dk = 65'=' (E (*) - (1 - *^) K (*)) (208) 

^*(l-*^)(l-^*^ + 12*^-^6-5^^K(*)d* = 


**'• 


( 3 * 2 - 1)55 
(l — *2) (l — 3*2 + *‘^) 
3*2 - 1 


K(fc) + E(*) 


-ifc4 


6 6 


/2 


- 2 6"5'= K(*)d* = * 


* ( 5*2 - 4*4 + *6 - 1 ) _1 


-1) 18 

(209) 

-E(*) 

(210) 


6 s'® 

*' 


,,, 6-5'=K(*)d* = 

((* 4-*2 + 1)K(*)-E(*)). 


( 211 ) 


5.3 Integrals from fragments of the equation for the com¬ 
plete elliptic integral of the second kind 

In equation ([5]) with y (*) a solution of equation (11551) we can take *(*) = ! 
and employ the recurrence relation: 


to obtain the integral 


dE(*) _ E(fc) -K(*) 
dk k 


J ^^rf*=K(*)-E(fc). 


( 212 ) 


(213) 


Using the identity (11861) equation (I155|) can be expressed in the alternative form: 


z" (*) + -/ (*)+ ! + 


1 - *2 


z (*) = 0 


and suitable fragmentary equations from equations (I155|) and (j214l) are: 


*"(*) + -*' (*) =0 
* 


1 


h" (fc) + (*) = 0 


(214) 

(215) 

(216) 
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h” (k) + yh' (k) + h = 0 
k 

h" {k)+h = Q 


(217) 

(218) 


\h'{k) + ^^h{k)=Q (219) 

+ (220) 

\h'{k)+h^Q. (221) 

k 

Equations (I215I) - (I221I) have the respective (non constant) solutions 


h (fc) = In (fc) 


( 222 ) 


hik) = k' 


Qcp-i {k)j 


h{k) = 


jjom 

\Yo{k)j 


h{k) = 


fsin(fc)'l 
\cos {k)j 


h{k) = k' 
h{k) = k'e^^^ 
hik)=e--^'^^ 


(223) 

(224) 

(225) 

(226) 

(227) 

(228) 


where in equation (12231) . ip = (l + -^/S) /2 is the Golden Ratio. The solutions 
(122211 - (I228|) give the respective integrals: 


fcln(fc)E(fc) 


dk = {l- In (fc)) E {k) + In (fc) K {k) 


(229) 


kk 


Ql-i (fc) 


E [k) dk = 


k'K (k) 


{pk^-l)E{k)\ (P^_,{k)\ {p-l)k 


k' 


Qi-i {k) 


k' 


-E(fc) 


P^^{k) 

Ql{k) 


J ^ Jo (fc) E (fc) dk = Jo (fc) (K (fc) - E (fc)) - kJi (k) E (k) 

E (k) dk = 

■“ l^CUS {K) j j 
\ C3in f 

^E(fc) - 


cos (fc) 1 k^ ( sin (k) 
— sin {k) 

cos (k) 


— sin (fc) 


k'"^ (cos {k) 

sin (fc) 
cos(fc) 


(E(fc)-K(fc)) 


(230) 

(231) 

(232) 


25 















fc'K (k) 


(233) 


/ 






E(fc) 

k' 



e5'='E(fc) dfc = 


(^fc'K(fc)- (^fc' + ^)E(fc)) 

(f) ) “ (^ + ^ 


(234) 

(235) 


5.4 Integrals involving the Golden Ratio 

Each transformed differential equation such as (I157p . (ll59p . (ll61l) and (11631) yields 
its own set of fragmentary equations and their corresponding integrals, similar 
in number to those givenabove for equations (11531) and (I155E and too numer¬ 
ous to present here. However, of particular interest are integrals related to 
equation (12301) above, which link Legendre/associated Legendre functions with 
E (fc), K {k) and the Golden Ratio (p. Equation (11631) above has the associated 
fragmentary equations: 


h" (k) - 


2k 


1-A:2 I l-fc2 (l_fc2)2 


h{k)=0 


h" {k) - 


2k 


1 - fc2 

which have the respective solutions: 


h' {k) 


1 


l-/c2 


h{k)=Q 


(236) 

(237) 


h (k) = 
h (k) = 


Pl-i (fc) 
Qi-i {k) 


’P^-iik) 

^Qtp—i {k) 

and substituting these results in equation ([2]) gives the respective integrals: 


(238) 

(239) 




ipk — 


3 — A:2 


2k 


E(A) + —K(A) 


P^-i(k) 

Q^-i (k) 




Qi (k) 


l-6P + k* (P^-iik) 


{kk')^^‘^ iQv-i {k) 


A\/k 

E(fc)dfc= 


k' 


(j)k — 


Z — k'" 
2k 


E(fc) + —K(fc) 

k 


Pv-i ik) 

^ip—\ (^) 




(240) 


(241) 
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6 Comments and conclusions 


A new Lagrangian method has been presented for deriving indefinite integrals 
of any function which obeys a second order linear ordinary differential equa¬ 
tion. The main result can also be proved very simply without variational calcu¬ 
lus. Various approaches have been presented to exploit the main formula and 
transformation methods have been given to multiply the number of interest¬ 
ing integrals obtainable. Sample results have been presented for some special 
functions, but these only scratch the surface of what is possible. The total 
number of indefinite integrals given in current tables and handbooks can be 
multiplied considerably with the method. All of the results presented here have 
been numerically checked using Mathematica [12]. 
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